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Abstract 

We consider magnetic fields generated by homogeneous isotropic and parity invariant turbulent flows. We show 
that simple scaling laws for dynamo threshold, magnetic energy and Ohmic dissipation can be obtained depending 
on the value of the magnetic Prandtl number. 
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Version francaise abregee II est a present admis que les champs magnetiques des etoiles voire meme des galaxies 
sont engendres par l'ecoulement de fluides conducteurs de l'electricite [1,2,3]. Ceux-ci impliquent des nombres de 
Reynolds cinetique, Re, et magnetique, R m , tres eleves (Re = VL/u, R m — poaVL, ou V est l'ecart-type des 
fluctuations de vitesse, L, l'echelle integrate de l'ecoulement, v, la viscosite cinematique du fluide, a, sa conductivite 
electrique et po, la permeabilite magnetique). Aucune experience de laboratoire ou simulation numerique directe 
des equations de la magnetohydrodynamique, ne permet 1' etude du probleme dans des regimes de parametres, 
Re et R m , d'interet astrophysique. II est done utile de considerer des hypotheses plausibles afin de pousser plus 
loin l'analyse dimensionnelle qui, a partir des parametres V, L, v, a, no et de la densite du fluide p, predit pour 
le seuil de l'effet dynamo et la densite moyenne d'energie magnetique, B 2 /2fio, saturee non lineairement au-dela 
du seuil, 

R c m = f(Re), (1) 

— = P V 2 g(R m ,Re). (2) 
Mo 

Dans le cas d'un ecoulement turbulent homogene isotrope, done de vitesse moyenne nulle, et invariant par 
symetrie plane, done sans helicite, les resultats des simulations numeriques les plus performantes realisees a ce 
jour montrent que R m augmente continuellement en fonction de Re [5]. Schekochihin et al. proposent que deux 
scenarios extremes, schematisees dans la figure 1, seront susceptibles d'etre observes lorsque les ordinateurs auront 
acquis la puissance requise pour effectuer des calculs a Re plus eleve : (i) une saturation R m — > constante, ou alors 
(ii) une croissance de la forme R m oc Re. D'autres simulations numeriques directes, realisees avec des ecoulements 
turbulents possedant un champ de vitesse moyen de geometrie fixee v(r) 7^ 0, semblent suivre le scenario (i) [13]. 
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Lorsque le nombre de Prandtl magnetique, P m — R m /Re — poav, est faible, P m <C 1, l'echelle de dissipation 
Joule du champ magnetique, l„ = LK^^ A , est grande par rapport a l'echelle de Kolmogorov Ik = LRe~ 3 / 4 . Le 
champ magnetique se developpe done a une echelle suffisamment grande pour ne pas etre affecte par la viscosite 
cinematique. Cette hypothese, couramment effectuee en turbulence, permet de conclure en faveur du scenario 

(i) . En effet, si v n'est pas pris en compte, l'analyse dimensionnelle impose R c m — > constante. II n'est done pas 
surprenant que les modelisations numeriques des grandes echelles, qui ne resolvent pas les echelles dissipatives, 
donnent ce resultat. Le scenario (i) sera done toujours observe a P m suffisamment faible sous reserve bien stir que 
Ton ait dynamo. 

II est cependant utile d'analyser le scenario (ii) d'autant plus que, comme nous pouvons le remarquer, il corres- 
pond a la prediction faite par Batchelor en 1950 [4] . En se basant sur une analogie entre l'equation de l'induction 
et celle de la vorticite, Batchelor avait estime que le seuil d'une dynamo engendree par un ecoulement turbulent 
devait correspondre a P m d'ordre unite, soit R^ oc Re. Meme si nous savons aujourd'hui que l'analyse de Bat- 
chelor est discutable, il est interessant de determiner sous quelle hypothese minimale sa prediction est correcte. 
Supposons done que nous nous limitions aux modes instables de champ magnetique, suffisamment localises au sein 
de l'ecoulement afin de ne pas etre affectes par les conditions aux limites. II est alors possible de ne pas prendre 
en compte l'echelle spatiale L, et l'analyse dimensionnelle impose pour le seuil, P m = constante, soit le scenario 

(ii) 7?^ oc Re. 

Les scenarios considered ci-dessus conduisent aussi a des predictions differentes pour la densite d'energie magne- 
tique engendree par effet dynamo. Le scenario (i) qui consiste a ne pas prendre en compte v revient a negliger 
la dependance en Re de g(R m ,Re) dans (2). Au voisinage du seuil, V est determine par R c m ~ p aVL et 
g(Rm) oc R m — R^ dans le cas d'une bifurcation supercritique. II en resulte [12] 

B 2 cx P (R m -R c m ). (3) 
pn(aL)^ 

Loin du seuil pour P m « 1, i?e » R m >• R c m , on peut supposer que B ne depend plus de a a condition que le 
champ magnetique se developpe a une echelle plus grande que l a . II en resulte alors l'equipartition entre energie 
magnetique et cinetique, B 2 / uo oc pV 2 , tel que suppose initialement par Biermann et Schliiter [15]. 

Un resultat completement different est obtenu dans le scenario (ii). II convient de considerer les parametres du 
probleme sous la forme equivalente, B, e = V s /L, L, v, a, uo et p. En effet, le champ magnetique a petite echelle 
est alimente par la puissance par unite de masse e qui cascade depuis l'echelle integrale, et il est done important 
de conserver ce parametre meme si Ton ne prend pas en compte explicitement L. L'analyse dimensionnelle conduit 
alors a 

— = pv^h(p m ) = h(p m ), (4) 

Mo VRe 



qui, pour P m ~ 1, n'est autre que le resultat obtenu par Batchelor en supposant que la saturation correspond a 
l'equipartition entre l'energie magnetique et l'energie cinetique a l'echelle de Kolmogorov. 

Revenons au cas P m <JC 1 qui correspond aux ecoulements de metaux liquides et plasmas a l'origine du champ 
magnetique des planetes et des etoiles (P m < 10 _ ). Dans ce cas, le champ magnetique se developpe a des echelles 
a priori comprises entre L et l a avec l a 3> Ik et il en resulte que R c m ne depend pas de variations de P m (ou de 
Re) et que B 2 / uo = pV 2 g(R m ) (scenario (i)). Interessons nous a la puissance dissipee par effet Joule par une telle 
dynamo. II faut a cet effet determiner a quelles echelles se developpe le champ magnetique. Utilisons pour cela un 
argument a la Kolmogorov en supposant que dans la zone inertielle, e'est a dire pour les nombres d'onde k tels 
que kl a <JC 1 <C kL, la puissance spectrale |B(fe)| 2 est independante de L, a and v. II en resulte 

\B\ 2 xpopeik'i. (5) 

Ceci n'est pas la seule possibility parmi les nombreuses predictions relatives au spectre de la turbulence magneto- 
hydrodynamique, mais dans le cas present, e'est probablement la plus simple. L'integration sur k redonne l'equi- 
partition B 2 /uo oc pV 2 . La contribution dominante a l'effet Joule provient de l'echelle l a - Nous obtenons 

— = - f \j\ 2 dk(x -4- [ k 2 \B\ 2 dk(x oc p^—, (6) 

a a J Mo " J Moo" L 

ou j est le vecteur densite de courant. Nous constatons done que la dissipation Joule est du meme ordre que 
la puissance totale disponible. Remarquons qu'il en serait de meme pour une dynamo de Batchelor suivant le 
scenario (ii) pour P m ~ 1, car bien que la densite d'energie soit plus faible, l'echelle caracteristique du champ 
magnetique Test egalement. 
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1. Introduction 

It is now believed that magnetic fields of stars and possibly galaxies are generated by the motion of electrically 
conducting fluids through the dynamo process [1,2,3]. These flows involve huge kinetic, Re, and magnetic, R m , 
Reynolds numbers (Re — VL/v, R m = pL^aVL, where V is the rms velocity amplitude, L is the integral length 
scale, v is the kinematic viscosity of the fluid, a is its electrical conductivity and po is the magnetic permeability). 
No laboratory experiments, neither direct numerical simulations are possible in the range of Re and R m involved 
in astrophysical flows. It is thus interesting to try to guess scaling laws for the magnetic field using some simple 
hypothesis. We consider here the minimum set of parameters, V, L, v, p®, a and p, the fluid density. We note that 
discarding global rotation makes our results certainly invalid for many astrophysical objects but not all of them. 
Rotation is indeed not assumed important for the galaxies which do not display a large scale coherent magnetic 
field [1,2,3]. Calling B its rms value, dimensional analysis gives 

R c m = f(Re), (7) 
for the dynamo threshold, and 
B 2 

— = pV 2 g(Rm,Re), (8) 
po 

for the mean magnetic energy density in the nonlinearly saturated regime. Our aim is to determine / and g in 
various regions of the parameter space (R m ,Re), assuming that turbulence is homogeneous, isotropic and parity 
invariant (thus with no mean flow and no mean magnetic field generation through an alpha effect). As already 
mentioned, this may look like an academic exercise compared to most natural dynamos. It is however not more 
academic that the concept of homogeneous and isotropic turbulence with respect to real turbulent flows. We thus 
expect that our simple arguments may shed some light on open problems concerning the effect of turbulence on 
the dynamo threshold and on the dynamic equilibrium between magnetic and kinetic energy. 

The dependence of the dynamo threshold R m — f(Re) in the limit of large Re is still an open problem, even 
in the case of a homogeneous isotropic and parity invariant turbulent flow. Note that parity invariance prevents 
the generation of a large scale magnetic field via an alpha effect type mechanism and isotropy implies zero 
mean flow. Recent direct numerical simulations show that R m keeps increasing with Re at the highest possible 
resolution without any indication of a possible saturation [5]. Schekochihin et al. thus propose that two limit 
scenarios, sketched in figure 1, could be observed when computers will be able to reach higher Re: (i) saturation, 
Rm ~ * constant, or (ii) increasing threshold in the form R m oc Re. 

A lot of work has been performed on the determination of R m as a function of Re for turbulent dynamos in the 
limit of large Re (or small P m ). We recall that (ii) has been proposed by Batchelor in one of the first papers on 
turbulent dynamos [4]. A lot of analytical studies have been also performed, mostly following Kazantsev's model 
[6] in order to show that purely turbulent flows can generate a magnetic field. Kazantsev considered a random 
homegeneous and isotropic velocity field, 5-correlated in time and with a wave number spectrum of the form k~ p . 
He showed that for p large enough, generation of a homogeneous isotropic magnetic field with zero mean value, 
takes place. This is a nice model but its validity is questionable for realistic turbulent flows. However, Kazantsev's 
model has been extrapolated to large Re. Various predictions, R m oc Re [7], R m — > constant « 400 for velocity 
spectra with 3/2 < p < 3 and no dynamo otherwise [8], or dynamo for all possible slopes of the velocity spectrum 
1 < p < 3 [9] have been found. These discrepancies show that extrapolation of Kazantsev's model to realistic 
turbulence cannot be rigorous. The calculation is possible only in the case of a 5-correlated velocity field in time, 
and 8(t — t'), which has the dimension of the inverse of time, should then be replaced by a finite eddy turn-over 
time in order to describe large Re effects. As already noticed, its choice is crucial to determine the behavior of 
R m versus Re. 

A different problem about turbulent dynamos has been considered more recently. It concerns the effect of 
turbulent fluctuations on a dynamo generated by a mean flow. The problem is to estimate to which extent the 
dynamo threshold computed as if the mean flow were acting alone, is shifted by turbulent fluctuations. This 
question has been addressed only recently [10] and should not be confused with dynamo generated by random 
flows with zero mean. It has been shown that weak turbulent fluctuations do not shift the dynamo threshold of 
the mean flow at first order. In addition, in the case of small scale fluctuations, there is no shift at second order 
either, if the fluctuations have no helicity. This explains why the observed dynamo threshold in Karlsruhe and 
Riga experiments [11] has been found in good agreement with the one computed as if the mean flow were acting 
alone, i.e. neglecting turbulent fluctuations. Recent direct numerical simulations have shown that in the presence 
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Figure 1. Dependence of the dynamo threshold R^ as a function of the Reynolds number Re. Scenario (i): R c m tends to a constant. 
Scenario (ii): R c m is proportional to Re. 

of a prescribed mean flow, v(r) =fc 0, R c m increases with Re at moderate Re but then seems to saturate at larger Re, 
thus following scenario (i). For the same flows, numerical modeling of large scales, large eddy simulations (LES) 
for instance, gives R^ ~ constant [13]. This last result follows from dimensional consideration as explained below, 
and has been also obtained for homogeneous isotropic turbulent non helical flows for which EDQNM closures have 
predicted R c m « 30 [14]. 



2. Turbulent dynamo threshold 

When the magnetic Prandtl number, P m = R m /Re — fioav, is small, P m <C 1, the Ohmic dissipative scale, 
l a — LR^^ 4 is much larger than the Kolmogorov Ik = LRe~ :i ^ 4 . Thus, if there is dynamo action, the magnetic 
field grows at scales much larger than Ik and does not depend on kinematic viscosity. This hypothesis is currently 
made for large scale quantities in turbulence and if correct, scenario (i) should be followed. If v is discarded, 
i?„ = constant indeed follows from dimensional analysis. It is thus not surprising that numerical models that do 
not resolve viscous scales, all gives this result, although the value of the constant seems to be strongly dependent 
on the flow geometry and on the model. We conclude that if dynamo action is observed for P m <C 1, the dynamo 
threshold is 

Rm — > constant when Re — > oo. (9) 

However, we emphasize that no clear-cut demonstration of dynamo action by homogeneous isotropic and parity 
invariant turbulence exists for P m <C 1. Experimental demonstrations as well as direct numerical simulations all 
involve a mean flow and analytical methods extrapolated to P m <C 1 are questionable. 

It may be instructive at this stage to recall the study on turbulent dynamos made more than half a century 
ago by Batchelor [4]. Using a questionable analogy between the induction and the vorticity equations, he claimed 
that the dynamo threshold corresponds to P m — 1, i.e. R c m oc Re, using our choice of dimensionless parameters 
(scenario (ii)). 

It is now often claimed that Batchelor's criterion P m > 1 for the growth of magnetic energy in turbulent 
flows is incorrect. However, the weaker criterion P m > constant (scenario (ii)) has not yet been invalidated 
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by direct numerical simulations or by an experimental demonstration without mean flow. It is thus of interest to 
determine the minimal hypothesis for which Batchelor's predictions for dynamo onset is obtained using dimensional 
arguments. To wit, assume that the dynamo eigenmodes develop at small scales such that the threshold does not 
depend on the integral scale L. Then, discarding L in our set of parameters, dimensional analysis gives at once 
P m — P^ — constant for the dynamo threshold, i. e. 

R c m oc Re. (10) 

It has been sometimes claimed that a non zero mean flow is necessary to get a dynamo following scenario (i). 
However, we note that even for a slow dynamo, i.e., growing on a diffusive time scale, the largest scales look 
stationary for a dynamo mode at wave length l a . For Kolmogorov turbulence, we indeed have, po&l 2 /(L/V) oc 
Rm Cl. This remains true for a k p spectrum for p < 3. 



3. Mean magnetic energy density 

Dimensional arguments can be also used to determine scaling laws for the mean magnetic energy density. For 
Pm <C 1 (scenario (i)), discarding v gives 

— = P V 2 go(R m ), (11) 
po 

where go is an arbitrary function. Close to threshold, the rms velocity V is given by poaVL ~ R m . In the case of 
a supercritical bifurcation, go (Rm) oc R m — R m , and we obtain [12] 

B 2 OC (Rm ~ Rm)- (12) 

Far from threshold, Re 3> R m ^> R m , ° ne could assume that B no longer depends on a provided that the magnetic 
field mostly grows at scales larger than l a . We then obtain equipartition between magnetic and kinetic energy 
densities, 

B 2 /p oc pV 2 , (13) 

as assumed by Biermann and Schliiter [15]. 

A completely different result is obtained in scenario (ii). Let us first recall that according to Batchelor's analogy 
between magnetic field and vorticity [4], the magnetic field should be generated mostly at the Kolmogorov scale, 
Ik = LRe~ 3/l4 , where the velocity gradients are the strongest. He then assumed that saturation of the magnetic 
field takes place for (B 2 ) / po oc pv K — pV 2 /V Re, where vk is the velocity increment at the Kolmogorov scale, 
t)| = y/ve. e = V 3 /L is the power per unit mass, cascading from L to Ik in the Kolmogorov description of 
turbulence. 

This can be easily understood, e = V 3 /L being the power per unit mass available to feed the dynamo, it may 
be a wise choice to keep it, instead of V in our set of parameters, thus becoming B, p, e, L, v, po and a. Then, if 
we consider dynamo modes that do not depend on L, we obtain at once 

— = pV^HPm) = h(P m ) (14) 



Mo VRe 

for saturation, where h(P m ) is an arbitrary function of P m - Close to dynamo threshold, P m ~ P m , we have 
h(Pm) oc Pm — P m if the bifurcation is supercritical. Only the prefactor pV 2 l\[~Re of (14) is the kinetic energy 
at Kolmogorov scale, that was assumed to be in equipartition with magnetic energy in Batchelor's prediction. 
This class of dynamos being small scale ones, it is not surprising that the inertial range of turbulence screens the 
magnetic field from the influence of integral size, thus L can be forgotten. We emphasize that a necessary condition 
for Batchelor's scenario is that the magnetic field can grow below the Kolmogorov scale, i.e. its dissipative length 
l a should be smaller than Ik, thus P m > 1. 

There is obviously a strong discrepancy between (13) and (14). The prefactors in these two laws are the upper 
and lower limits of a continuous family of scalings that are obtained by balancing the magnetic energy with the 
kinetic energy at one particular length scale within the Kolmogorov spectrum. It is not known if one of them is 
selected by turbulent dynamos. 
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4. Ohmic losses 

Ohmic losses due to currents generated by dynamo action give a lower bound to the power required to feed a 
dynamo. In order to evaluate them, it is crucial to know at which scales the magnetic field grows. Assuming that 
a dynamo is generated in the case P m <C 1 (scenario (i)), we want to give a possible guess for the power spectrum 
\B\ of the magnetic field as a function of the wave number k and the parameters p, e, L, u, /xo and a. Far from 
threshold, Re ^> R m 3> R c m , the dissipative lengths are such that Ik <SC la <§C L. For k in the inertial range, i.e. 
kl a « 1« kL, we may use a Kolmogorov type argument and discard L, a and v. Then, only one dimensionless 
parameter is left, and not too surprisingly, we get 



This is only one possibility among many others proposed for MHD turbulent spectra within the inertial range, 
but it is the simplest. Integrating over k obviously gives the equipartition law (13) for the magnetic energy. It is 
now interesting to evaluate Ohmic dissipation. Its dominant part comes from the current density at scale Z CT . We 
have 



We thus find that Ohmic dissipation is proportional to the total available power which corresponds to some kind 
of optimum scaling law for Ohmic dissipation. Although, this does not give any indication that this regime is 
achieved, we note that the above scaling corresponds to the one found empirically from a set of numerical models 
[16]. Their approximate fit, (B 2 /no)/(j 2 jo) oc L/V, indeed results from equations (15, 16). 
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